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We present numerical solutions of the semi-classical Boltzmann-Vlasov equation for fermion 
particle-antiparticle production by strong electric fields in boost-invariant coordinates in (1-1-1) and 
(3-1-1) dimensional QED. We compare the Boltzmann-Vlasov results with those of recent quantum 
field theory calculations and find good agreement. We conclude that extending the Boltzmann- 
Vlasov approach to the case of QCD should allow us to do a thorough investigation of how back- 
reaction affects recent results on the dependence of the transverse momentum distribution of quarks 
and anti-quarks on a second Casimir invariant of color SU(3). 



PACS numbers: 25.75.-q, 52.65.Ff, 12.38.Mh 

I. INTRODUCTION 

In recent papers, we have presented numerical quan- 
tum field theory calculations of the dynamics of fermion 
pair production by strong electric fields with back- 
reaction in boost-invariant coordinates in (1+1) and 
(3-1-1) dimensions [1, 2]. The purpose of the present pa- 
per is to compare these calculations with the results of 
numerical calculations using a semi-classical Boltzmann- 
Vlasov (BV) equation with a Schwinger source term for 
particle pair creation. We find that in (3-1-1) dimensions 
this semi-classical transport approximation works even 
better than it did in (1+1) dimensions. With the confi- 
dence that this model is working well for (3+1) dimen- 
sional quantum electrodynamics (QED), our program is 
to extend this calculation to quantum chromodynam- 
ics (QCD), where recently it has been shown that the 
WKB source term used by previous studies of pair pro- 
duction using the BV equation neglected an important 
term which depends on the second Casimir invariant of 
SU(3) [3, 4]. The BV equation is much quicker to imple- 
ment than the full field theory calculation and will let us 
explore the parameter space quickly before we perform 
more computer-intensive field theory calculations. 

The model we are using for the production of the par- 
ticles following a heavy ion collision is the so-called color 
flux tube model. The color flux tube model assumes that 
when two relativistic heavy ions collide multiple gluons 
are exchanged which leads to the formation of a strong 
color electric field. This model was studied extensively in 
the 1980's by several authors. These include Bialas, et al. 
[5, G, 7, H, 9, 10] and by Kajantie and Matsui [ ]. The 
idea of using a boost invariant Bolzmann-Vlasov equa- 
tion to study the time evolution of the plasma formed 



'Electronic address: john.dawson@unh.edu 
tElectronic address: bmihailaOlanl.gov 
fElectronic address: cooper@santafe.edu 



by the produced quarks and gluons was first put forward 
by Bialis and Czyz [">] and this was then generalized to 
include a Schwinger source term by Gatoff, Kerman and 
Matsui [i^]. At that time the validity of the BV approach 
was not known. However, once field theory calculations 
of this process were done in the 1990's [ ], it was clear 
that solving the BV equations with a Schwinger source 
term was a reasonable approximation. In the original 
work on QCD, the source term used was a WKB source 
term proposed by Casher, Neuberger and Nussinov [14], 
which recently has been shown to be incorrect by Nayak 
and collaborators [3, 4]. For constant chromoelectric 
fields the dependence on the second Casimir invariant 
can affect the transverse distribution of produced parti- 
cles by as much as 15% [1-j] which is a reason to correctly 
formulate the transport approach for the QCD plasma 
evolution and compare it to the field theory calculation. 

Our discussion of the BV equation in boost invariant 
coordinates for (3+1) dimensional QED follows closely 
in spirit work by Kluger, et al. [!(>, 17] and by Cooper, 
et al. [1 '] . We follow the method of solution used in these 
previous papers. In Section II, we discuss the classical 
theory for the boost-invariant coordinate system which 
we use in this paper and develop the equations needed 
for solutions of the BV equation. Numerical methods 
and results are discussed in Section III, and conclusions 
given in Section IV. 



II. CLASSICAL THEORY 

We wish to describe the dynamics of a relativistic par- 
ticle of mass M and charge e interacting with an elec- 
tromagnetic field in an arbitrary coordinate system. Let 
x'^^s) be the trajectory of a particle in space-time de- 
scribed parametrically by the arc- length ds, defined by 

{ds)^ = g^,,{x)dx^'dx^' . (2.1) 
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The velocity four- vector u'^(s) along the trajectory curve 
is given by 



(2.2) 



and the Lagrangian is 



1 



C[x^',u'']^ - M u^'{s)Uf,{s) + eu''{s)Af,{x) . (2.3) 



The canonical momenta Pp,{s) is given by 
dC 

^^^^^ ^ dvJ^ " + ' 



(2.4) 



where fcp(s) = Mm^(s) is the kinetic momentum. In 
terms of the kinetic momentum, Lagrange's equation give 



ds 



(2.5) 



where F^j^y{x) = d^Ay{x) — d„Af^{x) is the field tensor, 
which satisfies the Maxwell equations. 
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-g 



d^[^gF^'^{x)]^r{x) 



(2.6) 



where the current is the sum of convective and polariza- 
tion currents. The classical convective current is given 

by 

jconM(^)^g^ J dsu^'{s)S*[x~ x{s)] , (2.7) 

where x{s) is a solution of the equations of motion, and 
the sum goes over all species, particles, antiparticles, and 
spins. The energy momentum tensor densities for the 
particles t'^'^{x) and field Q^'^{x) are given by 

t^"'{x)^J2 J dsu''{s)k''{s)S*[x- x{s)] , (2.8a) 
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e^^'^ix) = - g^, F'-^'F^p + F^^g'^^Ffi, . (2.8b) 

The field energy-momentum tensor density satisfies 
1 



■.d^[^gQ''''{x)]=-F''"{x)J,{x). (2.9) 
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A. Trajectory solutions 

We next find trajectory solutions to the equations of 
motion in boost invariant coordinates. The Cartesian 
set of coordinates is designated by Roman letters: x°- = 
{t,x,y,z), with the metric rjab = diag( 1, — 1, — 1, — 1 ). 
Boost-invariant variables are designated by Greek letters: 
x^ = ( r, p, 0, 77 ) , where 



t = T cosh rj , z ~ T sinh r\ . 
X = p cos 6* , y = P sin 9 , 



(2.10) 



with the metric 

g^^(x) = diag( 1, -1, -p^, -r^ ) . 

The kinetic momentum in boost-invariant coordinates is 
then given by 



k-t^ — ( hP \ 
dx'^(s) 



(2.11) 



= M 



M{t',p',9',v'). 



Here a primed quantity means a derivative with respect 
to s. The mass shell restriction requires 

fc^fc^ = fc2 - fci - [fc,,/T]2 m2 , (2.12) 
where we have defined k± by 

kl = kl + [ko/pr = M^p'' + p'e'^]. 

So kr = LUk^,k^{T), where 

Wfc^,fc,(r) = ^Jkl + [k^/rY + M2 . (2.13) 

We restrict the vector potential and electric fields 
to be in the ry-direction and depend only on r, i.e.: 
A^{x) = (0,0,0,A^(r)). Then the only non- vanishing 
components of the field tensor are given by 

Fr.,^{x) = -F^^r{x) = drA^{T) = -tE{t) . (2.14) 

Here we have defined E{t) — — [3T-^,,(T)]/r. So then the 
Newton's Eqs. (2.5) become 



= eS(r)fc^(s)/r, 



M 



ds 
dknis) 



= eTE{T) kr{s) 



(2.15a) 
(2.15b) 



with kp and kg constants of the motion. Using the fact 
that kris) = Mdr/ds, Eq. (2.15b) becomes 



[A:„(T) + eA^(r)] = 0, 



(2.16) 



from which we conclude that p,, — kri{T) + ev4,,(r) is a 
constant of the motion. We can also define x- and y- 
components of the transverse momentum by 

k^ = kP cose - pk*^ sine (2.17) 

= M [p' cose- pO' sine] = Mx' , 

ky = kP sine + pk" cose (2.18) 

= M [p' sine + pe' cose] = My' . 

In cylindrical coordinates. 



k± cos ( 



k± sin ( 



(2.19) 



which defines the angle 0. By computing the Jacobians 
of these transformations, we show that volume elements 
are related by 



dkri-dky = k±dk± dcf) = 



dkp dkg 



(2.20) 
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B. Rapidity variables 

It will be useful to define rapidity momentum variables 
{r,y). These variables are defined by 



kt ^ r cosh y 



kz = r sinh y . 



(2.21) 



which can be related to our boost- invariant set (kr,kji) 

by 



k^ = M—^ 
as 



dXa 



from which we find 



kr = r cosh(?7 — y) 
k^/r = r sinh(?7 — y) . 



On the energy shell, we have 



ki 



'krj/r 



2 2 

~ r 



(2.22) 



(2.23a) 
(2.23b) 



(2.24) 



so on the energy shell, r = Mj_ = \/k'^ + . The 
Jacobian for this transformation is given by 



dk^/dr , dkr/dy 
dk^/dr , dkrj/dy 



so 



(2.25) 



(2.26) 



dk^ ( dkjj/r ) = rdrdy . 
We will use this result in Section II C below. 

C. The Boltzmann-Vlasov equation 



We define a particle distribution function f{x, k) such 
that the particle current density is given by (see, for ex- 
ample, Calzetta and Hu [ I ^]) 



Nf'ix) = / Bkk'' f{x,k) , 



(2.27) 



and the particle energy-momentum density tensor is 
given by 



f'^ix) = / mk^k" !{x,k) 



(2.28) 



where 



2RQ{ko)5{e-M^)d^k 
Dfc= /o N3 / — : 2.29 

with R a degeneracy factor. For a single species of 
fermions in (3-|-l)dimensions, counting particles, antipar- 
ticles, and spin, i? = 4. In a general coordinate system, 
the BV equation is given by (see for example Cooper, 
et al. [I.H] or Gatoff, et al. [12].) 



d 



where C{x, k) is a source term. Multiplying (2.30) by 
y/—g and integrating over Dfc gives 



1 



d^[V^N'^{x)]=C{x), (2.31) 



where 



C{x) = / Dfc fc° C{x, k) 



(2.32) 



So if C{x, k) — 0, particle number is conserved. Multi- 
plying (2.30) by fc" ^/—g and integrating over Dfc gives 



1 
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where 



dA^gt^^''{x)]-F-"{x)J„{x) = C''{x), (2.33) 



C'ix) = / Dfcfc'^fc'^C(a;,fc) , 



(2.34) 



and Ja{x) — eNa{x). So if C(a;,fc) — 0, combining 
Eqs. (2.9) and (2.33), we see that with no source term, 
the total energy-momentum tensor density. 



T^"'{x)^t^"'{x) + Q^"'{x) 



(2.35) 



satisfies a conservation law, T^'^;^(x) = 0. For our case, 
the source of particles is creation of particle-hole pairs 
via the Schwinger mechanism, so the particle number is 
not conserved and the energy-momentum tensor, using 
only convective currents, is also not conserved. 

In boost-invariant coordinates, we assume that the dis- 
tribution function is a function of f{T,k±,kjj) only. So 
choosing a surface element in the direction of constant 
T, we have dS = Td^x±dri where d'^x± = pdpdO is the 
perpendicular area, and 



R 



{2ny 



+00 

d^fc, / dfc, 



fc^'/(T,fc-^,fc-„) 

TcJfc^.feJr) 



(2.36) 

where d^fcj^ = fcj^dfc_L d0 and ujk±,k,,i''') is given by 
Eq. (2.13). The /i = component of (2.36) gives the 
number of particles per unit "volume" in boost-invariant 
coordinates: 



d^NjT) 

d'^x±dr] 



tNo{t) 



(2.37) 



R f f /'+00 



{2n) 

In terms of rapidity variables, Eq. (2.37) becomes 
R 



d^N 



{2^y 



[rd'^xx^drj] [d k±^dy] 



(2.38) 



X / '^^^^l^<5(r-M^)fc./(T,fc^,fc,) 



2r 



R 



{2nf 



[Td^x_Ld77] [d^k^dy]ujk^,k^ /(t, fc_L,fc,,) 
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So the momentum distribution in rapidity variables is 
given by 



fo(T,kJ 



d^A^ 



Rt 



Rt 

J2nf 



+00 
dfc„ 



drj 

dk„ 



O 

(2.39) 



But the only 77 dependence is through fc,,. Evaluating 
Eqs. (2.23) on the mass shell, we find 

Wfe^.fc^ cosh(77 - t/) , fc,, = r A/_l sinh(r/ - y) , 

(2.40) 

(2.41) 



so that for fixed y and k^, we find 
and (2.39) becomes 



R 



d'^x±d'^k±dy (2tt)^ 



+00 



dfc^/(r,fci,fc,), (2.42) 



and is independent of rapidity. 

Adding a Schwinger source term to the BV equation, in 
boost-invariant coordinates the only non- vanishing com- 
ponents of F^^{x) in our case are given in Eq. (2.14), so 
that Eq. (2.30) becomes 



d__ dA{T) d 
dr dr dk 



(2.43) 

where the source term is given by 

C(r, fc^, fc„) = P(r, k^) \eE{T)\ S{t, k^) S{kjT), (2.44) 
with P{T,kj_) a Pauli suppression factor evaluated at 

^T] — 0; 

P(T,fci) = l-2/(r,fci,0), (2.45) 
and S{t, kj_) is the Schwinger pair creation factor 



S{t, kj 



In 



^^{kl+M^}/\eE{T)\ 



(2.46) 



We solve Eq. (2.43) for f{T,kj_,kjj) using the method of 
characteristics. In Section II A, we found the particle tra- 
jectories and we showed that A:^(t) — prf — e yl^(r), where 
Pri is a constant of the motion. So the total derivative of 
/[r, A;_L, fc^(T)] with respect to t is given by 

d/[T, fc_L,fc,,(r)] 



df[T, k^,kr,{T)] _ JA^{t) 9/[t, k^,k,,{T)] 



dr 



dr 



dkn 



Assuming that f{TQ,k±,kri) — 0, we then have 

/(r,fc±,fc^)= f dr'r' P{T',k^)\eE{T')\S{T\k^) 

J To 

X S[kn + eA^ir) - eA^(r')] , (2.47) 
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FIG. 1: (Color online) Evolution of the distribution function 
/(r, fcxjO) for a typical case with M = 1, e = 1, A{to) = 
and E{to) = 4. 



which can be integrated to give 

/(r, k^,h,) = 51 [ 1 - 2/(t„, fci, 0) ] 5(r„, k^) (2.48) 

n 

X e(T„ -To)e(T-T„) . 

Here t„ are solutions of the equation 

k^ + e [ A(r) - A(t„) ] = , for To < t„ < r. (2.49) 

In order to step out f(T,k±,kri) as a function of t, we 
first solve (2.48) at k^ ~ 0, 

/(r,fci,0) = ^[l-2/(r„,fci,0)]5(T„,fci) (2.50) 

n 

X e(r„ - To) e(r - T„) , 

where now t„ is a solution of the equation A(t„) = A(t), 
for To < r„ < r. One such solution is for t„ = r. Se- 
lecting out this case, and setting 0(0) = 1/2, Eq. (2.50) 
becomes 



f{T,k^,0) = 



Sir, kj_)/2 + E,„<, P(t„, fci) ^(t„, fc^) 



l + S{T,k^) 

(2.51) 

We show a plot of /(r, fcj_, 0) in Fig. 1 for a typical case. 
With /(t, /cj^, 0) now known, we can solve Eq. (2.48) for 
the full /(r, fc_L, k^i). 



D. Maxwell's equations 

The only non- vanishing components of F^i, are given in 
Eq. (2.14), so Maxwell's equation (2.6) in boost invariant 
coordinates is given by 

drE{T) = tJ^{t) = -Mt)/t , (2.52) 
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with E{t) = —[drA{T)]/T. There are two types of cur- 
rents, convection currents arising from the flow of parti- 
cles and vacuum polarization currents, 



Mr) = J™"(t) -I- JP°Hr) 



(2.53) 



In a general frame, the convective current is given by the 
charge e times the ?7-component of the particle current 
N^' given in Eq. (2.27), 



J™"(t) = eN^ir) = ej Bkk^ /(r, fc^, fc,) . (2.54) 

Inserting the result for f(T,k±,kri) from Eq. (2.47) and 
integrating over /c^ gives 



eR 



X P{T',k^)\E{T')\S{T',k^) , (2.55) 



where we have put 

k,{T',T)=e[A{T')-A{T)] 



(2.56) 



The polarization current is determined by finding the ap- 
propriate current which, when added to the convection 
current, gives energy conservation. In Section IIP below, 
we found this current to be 



jr°\r)/r = sgn[EiT)] 



eR 
(2^ 



(2.57) 



X / A:_Ldfc_LWfc^,o(T)P(T, fc_L)S'(T,A:_L) 
Jo 

in agreement with Eq. (5.8) in Cooper et. al. [ ]. 



E. Particle creation 

The density of particles plus antiparticles at time r is 
given by Eq. (2.37). Substituting our solution (2.47) into 
this equation gives 



TiVo(T) = 



(2' 



— / k^dk^ dk^f{T,k^,k^) (2.58 

T^r Jo J-oo 

% f°ki_ dk^ / dr' r'P(T', k^) \E{t')\ Sir', k^) 

Jo Jrn 



(2' 



The rate of production of particles plus antiparticles can 
be obtained by differentiating (2.58) with respect to r 
and using the BV equation (2.43). This gives 

Tdr (27r)^ 

(2.59) 
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FIG. 2: (Color online) Proper-time evolution of the electro- 
magnetic fields A{t) and E{t) and the electric current J,j(r) 
for boost-invariant coordinates in (1-1-1) dimensions. Solutions 
of the BV equation are compared with results from the quan- 
tum field theory (QFT) calculation discussed in Ref. 1. Here 
we choose M = 1, e = 1, A(to) = and E{to) = 4. 



The particle production in terms of rapidity variables 
is obtained by substituting (2.47) into (2.42). This gives 



d^iV 



(2.60) 



d'^x±d'^k±dy 

= r' P{r', k^) \E{r')\ S{r', k^) . 

A picture of this distribution is shown in Fig. 8 as a 
function of r. Integrating (2.60) over kj_ gives 



TiVo(r) 



(2.61) 



1 diV 

17 d^ 

= / fc^dfc^/ dr'T'P{T',k^)\E{T')\S{r',kj^), 

(27r)^ Jo J TO 

where Aj_ is the perpendicular collision area. 

F. Energy-momentum tensor 

The field energy-momentum tensor density is given by 
(2.8b). For our case in boost- invariant coordinates, it is 
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T 

FIG. 3: (Color online) Proper-time evolution of the energy- 
momentum tensor (matter energy and longitudinal pressure) 
for boost-invariant coordinates in (1-1-1) dimensions. Parame- 
ters are the same as in Fig. 2 



diagonal and given by 




= ^ diag( E^,E^, p'^E^ ,~T^E^) 



(2.62) 
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T 

FIG. 4: (Color online) Time evolution of the electromagnetic 
fields A{t)/t and E{t) and the electric current J,)(r) for 
boost-invariant coordinates in (3-|-l)dimensions. Solutions of 
the BV equation are compared with results from the quantum 
field theory (QFT) calculation discussed in Ref. . Here we 
choose M = 1, e = 1, A(ro) = and E{to) = 4. 



The matter energy-momentum tensor is given by (2.28) 
which we write here as 



tpAr) = J Bkk^ k, f{T, fci, k^) (2.63) 
= dia.g{e,pp,p'^pe,T^Pjj). (2.64) 

So in the boost-invariant system, the matter energy and 
pressures are given by 



e = J Dkivl^j,^ f{T, k^,kr,) , 


(2.65a) 


Pp = j ^kklf{T, k^,kr,) , 


(2.65b) 


Pb = j Bk {ke/pY f{T, fc_L, fc^) , 


(2.65c) 


Pn = y Dfc (A:^/r)^ /(r, kj_,k,^) . 


(2.65d) 



Inserting the result for f{T,k±,kri) from Eq. (2.47) and 
integrating over k^ gives, for (2.65a) and (2.65d), 



eR /■+°° r /t'\ 

'^'■'^ (7) 

xP{T',k^)\E{T')\S{T',k^), (2.66a) 



eR 



+00 



\, Ak,{r',r)/r] 



xPiT',ki_)\E{T')\SiT',k^). (2.66b) 



where kri{T' ^t) and cjfej^(r',r) are given in Eqs. (2.56). 



Multiplying the BV equation (2.43) by a;| j, , 



and : 
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FIG. 5: (Color online) Time evolution of the matter energy 
(e), longitudinal pressure {p„), and transverse pressure (px) 
for boost-invariant coordinates in (3-|-l)dimensions, with the 
same parameters as in Fig. 4. 



tegrating over Dfc, gives 



/•oo /'i-oo 

■ / fcj^dfcj^ / dfc^ 

Jo J -oo 



+ eE{T) w/c_L,fe„(r) 



T dr 
dfjr, kj_,krj) 
dk„ 



(2.67) 



R \eE{T 



[ k^ dk^ ujk^.oir) P{r, S{t, 
T^r Jo 



For the first term in (2.67), we integrate by parts and 
note that 

T dr 



dr \ T 



So the first term becomes simply 
dr T 



(2.68) 



(2.69) 



For the second term in (2.67), we integrate by parts over 
kri and get 



E{t 



(r) [ 

— I Dkkr^ f{T,k±,kjj 



E{t) J™"(r) 



(2.70) 



where the convective current is given by (2.54). The last 
term in Eq. (2.67) can be written as 

E{r) .P,°\t)/t , (2.71) 
where the polarization current JP°'(r) is given by 

JP°'(r)/r = sgn[£;(r)]^ (2.72) 

X / A:_Ldfc_LWfc^,o(T)P(r, fc_L)S'(T,A;_L) , 



which is what we quoted in Eq. (2.57) in Section II D. 
Combining the results in Eqs. (2.69), (2.70), and (2.71), 
and noting that the total current is given by J-qiT) = 
J^°^^{t) + J^°^{t), the energy density and longitudinal 
pressure conservation equation reads 



de e+p, _E{t)J^{t) _ d[E'' /2] 



dr 



dr 



(2.73) 



where we have used Maxwell's equation (2.52). The total 
energy density and longitudinal pressure is given by 

£ = e + £;2/2 , V,^p,-E^/2. (2.74) 

Multiplying Eq. (2.73) by t gives an equation of state: 

dr{T£)+r, = . (2.75) 

The Pp{t) and pe{T) pressures are equal. We can prove 
this by introducing the integration variables kg — kg/ p 
and fc^ = krj/T, and putting 



Dfc = i? 
= R 
= R 
where now 



fc_L dfc_L /■+°° dkn 1 



27r 



(2.76) 

27r TWfc, ,k„ 



dfcp Akg dfc,, 1 



(27r)3 Tpuk^-k^^k^ 
dfcp dkg dkri 1 



'kl + kl + kl + M2 . (2.77) 
From (2.65b) and (2.65c), we find for the pressures 
'•+°° dfcp dfce dfc,, klf{T,kp,kg,kn) 



Pp{t) = R 



(27r)3 
(27r)3 



' —oo 

Pg{T) = Rlll ,^.3 — 
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FIG. 6: (Color online) Time evolution of the particle plus 
antiparticle density distribution function for boost-invariant 
coordinates in (1+1) dimensions, AN{T)/Ay, and in (3+1) 
dimensions, {l/As_)AN{T)/&y, respectively. In (3+1) dimen- 
sions, the BV calculation predicts slightly more particle pro- 
duction. 

SO Pp(t) = P0{t), as we claimed. Including the field pres- 
sure, we see that the total pressures also satisfy the re- 
lation Vp{t) — Vb^t), as required by conservation of the 
energy-pressure tensor. 

For the transverse pressure, we have 

p^{T)^Pp{T)+pe{T)^ Jmklf{T,k^,kr,) (2.78) 

XPiT',k^)\EiT')\SiT',k^). 

The shear pressure vanishes. In the next section, we com- 
pare results of solving the BV equation with a quantum 
field theory calculation in both (1+1) and (3+1) dimen- 
sions. 



III. NUMERICAL METHODS AND RESULTS 

The numerical procedure is as follows: we set i? = 2 
and i? = 4 for (1+1) and (3+1) dimensions, respectively, 
and choose units such that h = 1. We set M = 1 and 
e = 1, and take tq = 1/M = 1. Then we set A{tq) — 0, 
and choose a value of Eq. We set up a grid of values 
of between and fcfc^max, and compute initial values 
for Jjj{tq) and /(To,fcj_,0). We also compute a value for 
dJrj/dr at tq. We can then take a fourth-order Runge- 
Kutta step in r to find new values for A{t) and E{t), 
using a linear interpolator for values of Jri^T). We then 
compute values for f{T,kj_,0) from Eq. (2.50), J,,(r) us- 
ing Eq. (2.55), and dJ^/dr at the new value of r, and 




K 

FIG. 7: (Color online) Long proper-time (r — 200) 
transverse particle plus antiparticle distribution function, 
{l/A±)d^N{T = 200)/d^fcxdy, for boost-invariant coordi- 
nates in (3+1) dimensions. 



(1/A_^) d3N(T)/(d2k^dy) 




T 



FIG. 8: (Color online) Time evolution of the transverse par- 
ticle plus antiparticle distribution function, N{T)/d'^k±dy, 
for boost-invariant coordinates in (3+1) dimensions. See also 
Ref. . 



proceed in this way until we reach the final value of r. 
This method does not require computation of the full 
function /(r, fcj^,A:^) at the expense of an additional in- 
tegral over r; however /(t, fcj^,fc^) can be computed at 
any point along the way. 

We consider the case when Eq = 4, and compare the 
Boltzman-Vlasov (BV) results with two sets of recent 
quantum field theory (QFT) calculations done by us in 
(1+1) and (3+1) dimensional QED [ , -']. Comparisons 
for (l+l)-dimensional QED are shown in Fig. 2 for A{t), 
E{t), and J(t), and in Fig. 3 for components of the 
energy-momentum tensor. The BV calculation misses 
the fine structure noticed in the oscillations of the QFT 
electric current calculation, which has some features of 
quantum tunneling in a two- well potential (see Ref. 20), 
but otherwise is close in magnitude. The two calculations 
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get out of phase for large times, but this does not affect 
the calculation of the particle production which is dom- 
inated by the early-time dynamics. The BV calculation 
predicts larger energy density and longitudinal pressure, 
but about the same ratios of energy density to pressure. 

In Figs. 4 and 5 we compare the BV results to the 
QFT results for (3-|-l)-dimensional QED. We note that 
the fields, currents, energy density, and pressures all track 
very well together. The agreement between the BV and 
QFT calculations is better in (3-1-1) than in (1-fl) di- 
mensions, suggesting that the extra degrees of freedom 
perform some smoothing. In (3-1-1) dimensions we do not 
observe dephasing between the BV and QFT results at 
late times, at least as far as our calculations were carried 
out. Again, the BV calculation predicts larger values 
for the energy density and longitudinal pressure, but the 
transverse pressure for both calculations are fairly close 
to each other. There is no fine structure present in the 
(3-1-1) QFT results for the electric current, as discussed 
in Ref. 2. 

In Fig. 6, we show the particle plus antiparticle produc- 
tion per unit rapidity for the two calculations. In (1-1-1) 
dimensions, the particle plus antiparticle production per 
unit rapidity for both calculations are very close, aside 
from the fine structure. In (3-1-1) dimensions, the BV cal- 
culation predicts a slightly larger production than in the 
QFT results, which is consistent with the fact that the 
BV electric current depicted in Fig. 4 is slightly larger 
than the QFT current. Just, as in QFT, particles are be- 
ing created corresponding to the field gradients, with the 
major contribution coming from the initial field gradient. 
Subsequent smaller step increases are observed before the 
particle density saturates. 

Comparison of the late time (r = 200) transverse par- 
ticle plus antiparticle distributions for the BV and QFT 
calculations are shown in Fig. 7. The results are very 
close. Finally, in Fig. 8, we show the BV calculation for 
the entire time evolution of the transverse particle plus 
antiparticle distribution, which very similar to the one 
reported in Ref. 2, except for an approximate 5-10% dif- 
ference in magnitude. For r greater than about 80, there 
is no appreciable change in the shape of the distribution 
function, as expected, since by that time all particles have 
been produced by the field. 

IV. CONCLUSIONS 

We have presented here results of a non-equilibrium 
BV calculation of the time evolution of the quasiparti- 
cle distribution function for quarks in the presence of a 
proper-time evolving electric field with a Schwinger pair 



creation term in boost-invariant coordinates in (1-1-1) and 
(3-1-1) dimensions. We have then compared these re- 
sults with recent QFT calculations. Our one-dimensional 
results agree with previous results in Ref. and give 
reasonable agreement with the field theory calculations 
when short time scales are averaged over. What is ini- 
tially surprising is that in (3-1-1) dimensions, the short 
time scale fluctuations are not present in the field theory 
calculations so that agreement between the exact and the 
BV approximation for many macroscopic variables such 
as the time evolution of the electric field and the effec- 
tive energy density and pressures are quite good. The two 
methods differ in the particle production rate by about 
5-10%, mostly at low momentum transfers. 

It is at first quite surprising that the BV results are 
so close to the QFT results. A first-principles approach 
to deriving a BV-like equation for the exact field the- 
ory equations in scalar electrodynamics in (1-1-1) dimen- 
sion has been given in Ref. where obtaining a local 
Vlasov source term from the non-local equation for the 
adiabatic number operator seemed to follow from phase 
decoherence of the quantum density matrix. In (3-1-1) 
dimensions we would imagine that this phase decoher- 
ence would occur more quickly than in (1-1-1) dimensions, 
which would make the quantum to classical transition 
quite rapid. This would then be the reason why the semi- 
classical approach presented here works better in (3-1-1) 
than in (H-1) dimensions. 

The fact that the BV calculations are computation- 
ally much faster than solving the field theory equations 
makes them a good candidate for extending this work 
to the case of QCD, where the computer time required 
for a full QFT calculation can become prohibitive for 
an exhaustive investigation of the two SU(3) Casimir in- 
variants parameter space. If the BV approach with the 
correct Schwinger source term proves to be as accurate in 
QCD as in QED then it would be very helpful in explor- 
ing parameter space so that the Casimirs dependence for 
the transverse distribution function can be better under- 
stood for the case when back-reaction is included. We 
intend to explore this possibility in a subsequent publi- 
cation. 
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